An explicit harmonic code for black-hole evolution using excision

We describe an explicit in time, finite-difference code designed to simulate black holes by using the excision method. The code is based upon the harmonic formulation of the Einstein equations and incorporates several features regarding the well-posedness and numerical stability of the initial-boundary problem for the quasilinear wave equation. After a discussion of the equations solved and of the techniques employed, we present a series of testbeds carried out to validate the code. Such tests range from the evolution of isolated black holes to the head-on collision of two black holes and then to a binary black hole inspiral and merger. Besides assessing the accuracy of the code, the inspiral and merger test has revealed that individual apparent horizons can touch and even intersect. This novel feature in the dynamics of the marginally trapped surfaces is unexpected but consistent with theorems on the properties of apparent horizons

Testing numerical relativity with the shifted gauge wave

Computational methods are essential to provide waveforms from coalescing black holes, which are expected to produce strong signals for the gravitational wave observatories being developed. Although partial simulations of the coalescence have been reported, scientifically useful waveforms have so far not been delivered. The goal of the AppleswithApples (AwA) Alliance is to design, coordinate and document standardized code tests for comparing numerical relativity codes. The first round of AwA tests have now being completed and the results are being analyzed. These initial tests are based upon periodic boundary conditions designed to isolate performance of the main evolution code. Here we describe and carry out an additional test with periodic boundary conditions which deals with an essential feature of the black hole excision problem, namely a non-vanishing shift. The test is a shifted version of the existing AwA gauge wave test. We show how a shift introduces an exponentially growing instability which violates the constraints of a standard harmonic formulation of Einstein's equations. We analyze the Cauchy problem in a harmonic gauge and discuss particular options for suppressing instabilities in the gauge wave tests. We implement these techniques in a finite difference evolution algorithm and present test results. Although our application here is limited to a model problem, the techniques should benefit the simulation of black holes using harmonic evolution codes

The Initial-Boundary Value Problem in General Relativity

In this article we summarize what is known about the initial-boundary value problem for general relativity and discuss present problems related to it.Comment: 11 pages, 2 figures. Contribution to a special volume for Mario Castagnino's seventy fifth birthda

Lithologic influences on groundwater recharge through incised glacial till from profile to regional scales: Evidence from glaciated Eastern Nebraska

Variability in sediment hydraulic properties associated with landscape depositional and erosional features can influence groundwater recharge processes by affecting soil-water storage and transmission. This study considers recharge to aquifers underlying river-incised glaciated terrain where the distribution of clay-rich till is largely intact in upland locations but has been removed by alluvial erosion in stream valleys. In a stream-dissected glacial region in eastern Nebraska (Great Plains region of the United States), recharge estimates were developed for nested profile, aquifer, and regional scales using unsaturated zone profile measurements (matric potentials, Cl 2 and 3 H), groundwater tracers (CFC-12 and SF 6 ), and a remote sensing-assisted water balance model. Results show a consistent influence of till lithology on recharge rates across nested spatial scales despite substantial uncertainty in all recharge estimation methods, suggesting that minimal diffuse recharge occurs through upland glacial till lithology whereas diffuse recharge occurs in river valleys where till is locally absent. Diffuse recharge is estimated to account for a maximum of 61% of total recharge based on comparison of diffuse recharge estimated from the unsaturated zone (0–43 mm yr 21 ) and total recharge estimated from groundwater tracers (median 58 mm yr 21 ) and water balance modeling (median 56 mm yr 21 ). The results underscore the importance of lithologic controls on the distributions of both recharge rates and mechanisms

Exploiting gauge and constraint freedom in hyperbolic formulations of Einstein's equations

We present new many-parameter families of strongly and symmetric hyperbolic formulations of Einstein's equations that include quite general algebraic and live gauge conditions for the lapse. The first system that we present has 30 variables and incorporates an algebraic relationship between the lapse and the determinant of the three metric that generalizes the densitized lapse prescription. The second system has 34 variables and uses a family of live gauges that generalizes the Bona-Masso slicing conditions. These systems have free parameters even after imposing hyperbolicity and are expected to be useful in 3D numerical evolutions. We discuss under what conditions there are no superluminal characteristic speeds

Geometrical optics analysis of the short-time stability properties of the Einstein evolution equations

Many alternative formulations of Einstein's evolution have lately been examined, in an effort to discover one which yields slow growth of constraint-violating errors. In this paper, rather than directly search for well-behaved formulations, we instead develop analytic tools to discover which formulations are particularly ill-behaved. Specifically, we examine the growth of approximate (geometric-optics) solutions, studied only in the future domain of dependence of the initial data slice (e.g. we study transients). By evaluating the amplification of transients a given formulation will produce, we may therefore eliminate from consideration the most pathological formulations (e.g. those with numerically-unacceptable amplification). This technique has the potential to provide surprisingly tight constraints on the set of formulations one can safely apply. To illustrate the application of these techniques to practical examples, we apply our technique to the 2-parameter family of evolution equations proposed by Kidder, Scheel, and Teukolsky, focusing in particular on flat space (in Rindler coordinates) and Schwarzchild (in Painleve-Gullstrand coordinates).Comment: Submitted to Phys. Rev.

An objective method for determining principal time scales of coherent eddy structures using orthonormal wavelets

A new, parameter-free method, based on orthonormal wavelet expansions is proposed for calculating the principal time scale of coherent structures in atmospheric surface layer measurements. These organized events play an important role in the exchange of heat, mass, and momentum between the land and the atmosphere. This global technique decomposes the energy contribution at each scale into organized and random eddy motion. The method is demonstrated on vertical wind velocity measurements above bare and vegetated surfaces. It is found to give nearly identical results to a local thresholding approach developed for signal de-noising that assigns the wavelet coecients to organized and random motion. The eect of applying anti-and/or near-symmetrical wavelet basis functions is also investigated.

Semi-automatic segmentation of subcutaneous tumours from micro-computed tomography images

Cataloged from PDF version of article.This paper outlines the first attempt to segment the boundary of preclinical subcutaneous tumours, which are frequently used in cancer research, from micro-computed tomography (microCT) image data. MicroCT images provide low tissue contrast, and the tumour-to-muscle interface is hard to determine, however faint features exist which enable the boundary to be located. These are used as the basis of our semi-automatic segmentation algorithm. Local phase feature detection is used to highlight the faint boundary features, and a level set-based active contour is used to generate smooth contours that fit the sparse boundary features. The algorithm is validated against manually drawn contours and micro-positron emission tomography (microPET) images. When compared against manual expert segmentations, it was consistently able to segment at least 70% of the tumour region (n = 39) in both easy and difficult cases, and over a broad range of tumour volumes. When compared against tumour microPET data, it was able to capture over 80% of the functional microPET volume. Based on these results, we demonstrate the feasibility of subcutaneous tumour segmentation from microCT image data without the assistance of exogenous contrast agents. Our approach is a proof-of-concept that can be used as the foundation for further research, and to facilitate this, the code is open-source and available from www.setuvo.com. © 2013 Institute of Physics and Engineering in Medicine